O/H in the Local Bubble 



Cristina M. Oliveira 2 , Jean Dupuis 2 , Pierre Chayer 2,3 , and H. Warren Moos 2 



ABSTRACT 

We present new measurements of the oxygen gas-phase abundance 
along four sightlines probing gas inside the Local Bubble, using data ob- 
tained with the Far Ultraviolet Spectroscopic Explorer (FUSE) and the 
Hubble Space Telecope (HST). Combining these with seven previously pub- 
lished iV(OI) and iV(HI) measurements we determine a mean O/H ratio for the 
Local Bubble, (0/H) LB = (3.45 ± 0.19) xl(T 4 (la in the mean). Our result is 
virtually identical to O/H = (3.43 ± 0.15) xlO -4 derived from data for sightlines 
probing gas well outside the Local Bubble. In contrast to the D/O and D/H 
ratios, which seem to have different values beyond the Local Bubble, our results 
show that the O/H ratio is constant and single- valued both inside and outside the 
Local Bubble, for low density sightlines, with column densities up to log iV(H) ~ 
21. In addition, the similarity of the two values above suggests that the net effect 
of the uncertain 01 /-values in the FUSE bandpass is not significant. Taking 
into account the latest determinations of the Solar abundance of oxygen, our 
result implies that ~25% of the oxygen in the Local Bubble is depleted onto dust 
grains. The similarity of the value for O/H along low average density sightlines 
in the Local Bubble with that of denser sightlines beyond may permit a limit on 
the diluting effects of infalling low metallicity gas. 

Subject headings: ISM: Abundances — ISM: Evolution — Ultraviolet: ISM - 
Stars: Individual (WD 0050-332, WD 0232+035, WD 1254+223, WD 2152-548) 
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1. INTRODUCTION 

The abundance of oxygen in the interstellar medium (ISM) is an important parameter 
of galactic chemical evolution models. In addition, oxygen, the third most abundant element 
in the ISM after hydrogen and helium, plays an important role in interstellar chemistry and 
in dust formation models. 

? used the weak intersystem A1356 01 transition to determine 0/H along the line of 
sight to 13 stars (with distances ranging from 130 to 1500 pc), with log iV(HI) = 20.18 - 
21.28 and log /(H 2 ) = -5.21 - -0.20. They found a remarkably constant O/H = (3.43 ± 
0.15)xl0 -4 (where the / value has been adjusted; see §6). ? used the same OI transition 
to determine O/H toward 19 stars (with distances ranging from 0.84 to 5.01 Kpc) with log 
/(H 2 ) = -1.52 to -0.33 and logiV(HI) = 20.94 - 21.41, deriving O/H = (4.08 ± 0.13) xl0~ 4 
(la in the mean). More recently, ? used the same transition to determine O/H toward 36 
new sightlines, with distances ranging from ~0.5 to 6.5 Kpc. Combining the new results with 
previously published HST results they find a bimodal distribution in O/H, with O/H = (3.90 
± 0.11) xlO" 4 along low density sightlines (n H < 1 cm" 3 ) and O/H = (2.84 ± 0.12) xlO -4 
along denser sightlines. In addition, these authors find, for the low density sightlines (rin < 
1 cm -3 ), that the O/H ratios diverge when one considers pathlengths shorter or longer than 
800 pc. They find (0/H) short = (3.47 ± 0.07) xlO" 4 , vs, (0/H) long = (4.37 ± 0.09) xlO" 4 . 

In the Local Bubble, iV(OI) is low enough that the A1356 transition is undetected; 
its /-value is ~13 times weaker than that of the weakest, measurable, OI transition in 
the FUSE bandpass (A974). For sightlines within the Local Bubble A^(OI) can only be 
determined with FUSE data. The FUSE bandpass contains numerous O I transitions, with 
different oscillator strengths, which allow accurate OI column densities to be determined. 
The strong A1302 transition (/ = 4.80 xlO -2 ), available only with STIS, is always saturated, 
and as such, oxygen column densities derived solely from this transition are not reliable. 

The studies by ? and ? indicated that O I is a good tracer of H I in the nearby Galactic 
disk. However, only a few sightlines inside the Local Bubble have published O/H ratios (see 
??, and references therein). 

In this paper, we use data obtained with FUSE and the Space Telescope Imaging 
Spectrograph (STIS), onboard HST, in conjunction with previously published OI and HI 
column densities, to determine the gas-phase O/H ratio in the Local Bubble, (0/H) L b- The 
observations and data reduction are described in §2. Determinations of A^(OI) and A^(HI) 
are presented in §3 and §4, respectively. The results are presented in §5 and discussed in §6. 
Our finding are summarized in §7. All results are quoted at the la level, unless otherwise 
noted. 
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2. OBSERVATIONS AND DATA PROCESSING 

We use data obtained with FUSE to determine N(0 1) and data obtained with STIS on- 
board HST to determine the neutral hydrogen column densities. The stellar properties of 
the 4 stars are summarized in Table 1. The observations are described with details below. 



2.1. FUSE Observations 

The FUSE observatory consists of four co-aligned prime-focus telescopes and Rowland- 
circle spectrographs that produce spectra over the wavelength range 905 - 1187 A with a 
spectral resolution of ~ 15 - 20 km s -1 (wavelength dependent), for point sources. Two 
of the optical channels employ SiC coatings, providing reflectivity in the wavelength range 
~ 905 - 1000 A, while the other two have LiF coatings for maximum sensitivity above 
1000 A. Dispersed light is focused onto two photon-counting microchannel plate detectors. 
The FUSE mission, its planning, and on-orbit performance are discussed by ? and ? 

Table 2 summarizes the FUSE observations of the four targets studied in this work. 
The data were obtained through both the large (LWRS, 30" x 30") and medium sized aper- 
tures (MDRS, 4" x 20") in histogram and time-tagged modes (HIST and TTAG, respectively 
in Table 2). The two-dimensional FUSE spectra are reduced using the CalFUSE pipeline 
V2.4.1 1 . The processing includes data screening for low quality or unreliable data, ther- 
mal drift correction, geometric distortion correction, heliocentric velocity correction, dead 
time correction, wavelength calibration, detection and removal of event bursts, background 
subtraction, and astigmatism correction. The spectra are aligned by cross-correlating the 
individual exposures over a short wavelength range that contains prominent spectral features 
and then coadded by weighting each exposure by its exposure time, using the CORRCAL 
software developed by S. Friedman. All the spectra are binned to three pixel samples, or 
~20 mA, for analysis (the line spread function, LSF, is about 11 pixels or ~70 mA wide). 
For each target, observations obtained through the same aperture are coadded in order to 
increase the S/N of the dataset. 



1 The CalFUSE pipeline reference guide is available at http://fuse.pha.jhu.edu/analysis/pipeline_reference.html 
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2.2. STIS Observations 

Table 3 summarizes the STIS data for the stars. The data are reduced using the standard 
STSDAS pipeline within IRAF. The CALSTIS pipeline (version 2.3 ?) is used to reduce the 
data obtained with the E140H and E140M echelles. The echelle spectral orders in the regions 
adjacent to Lya are combined using a weighted averaging scheme where the orders overlap. 
For stars where there are multiple observations the different exposures are combined by 
weighting them by the exposure time, to increase the S/N of the data. We use a single- 
Gaussian with a FWHM of 1.2 and 1.3 pix to describe the line spread function of the E140H 
and E140M data, respectively. 



3. Determination of A^(0 I) 

We use apparent optical depth, curve of growth and profile fitting methods (hereafter 
AOD, COG, and PF, respectively) to derive N(Ol). We follow the procedures outlined in ? 
to use these techniques, with the following modification: a single- Gaussian with a FWHM of 
10.5 pixels (constant across all wavelengths and channels and during the fitting procedure) 
is used to describe the FUSE line spread function. Table 4 summarizes the wavelengths and 
/-values (?) of the 01 transitions used in this study, as well as their measured equivalent 
widths. Table 5 presents, for each star, N(0 1) derived with the different methods, as 
well as the adopted iV(OI) values. We indicate in parentheses which line was used to 
determine the AOD result which is quoted in Table 5. For the AOD measurements we use 
the weakest O I line detected along each sightline, when the derived column densities have 
reasonable uncertainties, such as for WD 0050-332 and WD 1254+223. For WD 0232+035 
and WD 2152— 548, the AOD results quoted are from stronger lines. For these sightlines, 
AOD measurements on the weakest lines are consistent with those of the stronger lines 
used, which provide smaller uncertainties in the column densities. The adopted N(0 1) 
values represent compromises between the O I column densities obtained with the different 
methods. We chose to, subjectively, adopt N(0 1) values that are the midpoint of the range of 
measurements obtained with the different methods, with uncertainties that overlap with the 
extreme measurements. We note that no molecular hydrogen was detected along any of the 
four sightlines presented here. Below we comment briefly on each line of sight independently. 
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3.1. WD 0050-332 

The weighted average of the AOD measurements (on the SIB and S2A channels) on 
A924.9 yields log iV(OI) = 15.15 ± ° f 2 , consistent with AOD on A929.5 (log JV(OI) = 15.18 
t o°ol)- The COG fit yields log N(0 1) = 14.98 + g;{g. The STIS E140H data for this sightline 
shows at least two components in NI and Sill separated by ~10 km s -1 . The velocity of 
the stronger component, v & = 6.8 km s _1 , is consistent with the expected velocity of the 
LIC along this direction (?), i> L ic = 6.2 ± 1.2 km s -1 . Using the double-component velocity 
structure derived from these species, and fitting N I, Si II and O I together we find log N(0 1) 
= 15.09 ± 0.09. We adopt then log iV(OI) = 15.07 ± 0.09. We use more data than ? which 
find log iV(OI) = 15.26 + 



3.2. WD 0232+035 

AOD measurements on the A976 OI transition yields log jV(OI) = 15.01 ± 0.03. The 
COG fit to O I transitions only in the FUSE bandpass yields log iV(OI) = 15.06 t aoe- 
When the A1302 O I transition (present in the STIS data) is included in the COG fit we 
find log A(OI) = 14.96 t. q Q3- ? have analyzed the STIS data for this star. They find a 
stronger component at v Q +3.1 km s _1 , and a weaker component at v & +17.6 km s -1 . The 
latter component is identified with the LIC. We use the two-component velocity structure 
derived by these authors to fit the O I transitions in the FUSE bandpass. Using PF we derive 
log iV(OI) = 14.94 t o;o4- Figure 1 presents fits to some of the OI lines used with profile 
fitting, with the two-component model discussed above. Dotted lines represent the fit prior 
to convolution with the instrument LSF, and dashed lines indicate the adopted continua. 
We adopt then log iV(OI) = 15.00 ± 0.06 for this sightline, as a compromise between the 
AOD, COG and PF measurements. 



3.3. WD 1254+223 

The low column density along this line of sight, combined with the S/N of the data 
allows us to detect only the A988 (triplet) and A1039 OI transitions in the FUSE data. 
Since the A988 transition is strongly contaminated with terrestrial airglow, we use only the 
weak A1039 transition in conjunction with the strong A1302 line (from the STIS data) to 
constrain A^(OI) along this line of sight. 

We find log iV(0 1) = 14.22 ± 0.03, from a weighted average of the AOD measurements 
(on four channels) on the A1039 transition. The COG fit to the A1039 and A1302 transitions 
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Fig. 1.— Fits to some of the OI lines used in PF for the WD 0232+035 sightline, with 
the two-component model discussed in the text. Dotted lines represent the fit prior to 
convolution with the instrument LSF, dashed lines and dot-dash lines represent the two 
components. Solid lines represent the adopted continua as well as the fit after convolution 
with the instrument LSF. 
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yields log JV(OI) = 14.27 ± 0.04. We adopt then log iV(OI) = 14.25 t Our result is 
consistent, within the uncertainties, with log iV(OI) = 14.38 1 q ^ from ?. Their result is 
based on the possible detection of the A976 transition which is ~3 times weaker than A 1039 
{ T o,u)39 — 0-24), not detected in our dataset. The difference between the two datasets might 
be due to the CalFUSE pipeline used; ? used CalFUSE v2.0.5, while we use v2.4.1. 



3.4. WD 2152-548 

We derive log N(Ql) = 15.35 1 ^H) using the AOD method on the A936 transition. 
The COG and PF techniques yield log iV(OT) = 15.52 ± and log JV(OI) = 15.30 + 
respectively. We adopt log iV(0 1) = 15.45 ± 0.10 in good agreement with log iV(0 1) = 15.49 
1 ° f 7 derived by ?. 



4. Determination of A^(H I) 

We use observations of the H I Lya transition obtained by the Space Telescope Imaging 
Spectrograph onboard HST in conjunction with stellar models, to determine A r (HI) along 
three of the four sightlines studied in this work. Below we discuss each star separately. The 
adopted A^(HI) are presented in Table 6. 



4.1. WD 0050-332 

4.1.1. Stellar Model 

Several authors have used the Balmer line series to determine T c r and log g for this 
white dwarf. Using pure-hydrogen LTE stellar models ? find T cff = 35,816 ± 148 K and 
log g = 7.679 ± 0.031, while ? find T cff = 36,320 ± 200 K and log g = 7.92 ± 0.03. 
? have also analyzed the HI Balmer lines, using non-LTE pure-hydrogen stellar models. 
These authors find T cS = 35,660 ± 135 K and log g = 7.93 ± 0.03. To take into ac- 
count the stellar H I absorption a set of nine pure-H NLTE stellar models corresponding to 
T cS = (35,000;36,000;37,000) K and log g = (7.65;7.80;7.95) was constructed. The range in 
the stellar parameters, T e s and log g, was chosen to include the published T e g and log g values 
in the literature, as well as their associated uncertainties. 

We use the velocity shift between photospheric and interstellar absorption determined by 
? to align the stellar model with the data. From an analysis of STIS E140M and IUE data, 
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these authors find Vpn — Vism ~ +25 km s x , corresponding to a shift of ~18 pix. 

4-1.2. Ly a Profile Fitting 

We follow the procedure outlined in §2.2 to calibrate and coadd the STIS data covering 
the Lya region. Profile fitting of the HI Lya transition alone without constraining the 
Doppler parameter, b, leads to unreasonable estimates of iV and b. We use then the Lyman 
edge, covered by FUSE in the SiClB channel, to constrain b which is then used as a constraint 
when fitting the Lya region. By fitting the HI Lyman series transitions between 913.0 A and 
919.2 A with the continuum fixed, we determine &hi = 8.0 1 ^ km s -1 . In all our fits we use 
the fitting code Owens (for a more detailed description of this code see ??). After aligning 
and normalizing the data by each of the stellar models described above fit the Lya region. 
No DI is included in the fit and the b- value is fixed at 8.0 km s -1 . In addition, a first-order 
polynomial is used to model the continuum after normalization by the model atmosphere, 
to take into account uncertainties in the instrument sensitivity. Figure 2 shows the best 
fit model, corresponding to T eff = 36,000 K and log g = 7.80, overplotted on the coadded 
data. We show also the two models which produce the most extreme H I column densities, 
which lead to our la systematic uncertainties. These are the models with T c g = 35,000 K, 
log g = 7.65 and T cff = 37,000K, log g = 7.95. We find log iV(HI)= 18.57 ± 0.02 ± where 
the first uncertainties quoted take into account the statistical uncertainties associated with 
using the best fit model (T e g = 36,000 K, log g = 7.80) and the second quoted uncertainties 
reflect the systematic uncertainties associated with using different stellar models. Using a 
conservative approach to combine the two uncertainties we adopt log iV(HI)= 18.57 t o!o4- 
The value of bm derived above, by fitting the HI Lyman series transitions, implies T max = 
3840 K, considerably lower than the mean temperature, T = 6680 ± 1490 K, derived by ? 
from an analysis of sightlines within 100 pc. To ascertain the influence of &hi on iV(HI) we 
performed multiple fits using the best fit stellar model and fixing bai to different values. We 
find that iV(HI) agrees at the la level with the value quoted above, for 6hi < 20 km s _1 . 
Considering typical b- values for HI along similar sightlines (?), 20 km s _1 is a reasonable 
upper limit for 6hi- Hence our derived iV(HI) should not be biased by the &hi adopted. 

Our adopted iV(HI) is in good agreement with the values determined by ? and ? 
(log iV(HI) = 18.49 ± ° f 2 and log iV(HI) = 18.51 ± ° f 3 , respectively) from analyses of 
EUVE data for this star. 
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Fig. 2. — HI Lya absorption for WD 0050— 332. The histogram shows the coadded and 
merged STIS E140H observations of this line of sight; the thick solid line shows the adopted 
stellar model, corresponding to T e s = 36,000 K and log g = 7.80. Also displayed are the two 
models which produce the lowest iV(HI) (T e g = 35,000 K, log g = 7.65, dashed line) and 
the highest iV(HI) (T cS = 37,000 K, log g = 7.95, dashed-dotted line). Geocoronal emission 
is annotated with ©. 
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4.2. WD 1254+223 

4.2.1. Stellar Model 

Determinations of the effective temperature and gravity for this white dwarf range from 
T eff = 38,686 K, log g = 7.662 (?) to T cS = 41,200 K, log g = 7.85 (?). 

To determine the influence of the different stellar models in the measured interstellar 
H 1 column density we constructed a grid of nine pure-H NLTE stellar models corresponding 
to the values T cff = (39,000; 40,000; 41,000) K and log g = (7.60; 7.75; 7.90). ? tentatively 
determined a velocity shift of vpu — vi SM = +21 km s _1 between the photospheric and the 
interstellar absorption. We use the shift determined by these authors to align the stellar 
models with the data. 



4-2.2. Lya Profile Fitting 

For each of the nine computed stellar models we normalize the data by the stellar model 
prior to profile fitting. The D I Lya transition is also clearly visible in the data (see Figure 3) 
so we fit both D I and H I. Because H I can sometimes trace gas not detected in other species 
we fit the Lya D I and H I profiles in separate absorption components. We use a second- 
order polynomial to model the continuum after normalization by the model atmosphere. 
The Doppler parameter is a free parameter of the fit and we derive bm — 11.6 1 q 5 km s_1 
and 6di — 7.6 1 J'y km s _1 . The 6-values derived here for HI and DI are consistent with the 
temperature and turbulent velocity (T = 7000 ± 2800 K, £ = 1.2 ± \\ km s" 1 ) derived for this 
line of sight by ?. Figure 3 displays the data overplotted with the best fit model (T cS = 40,000 
K and log g = 7.75) and the two models that produce the lowest iV(H I) (T cS = 39,000 K and 
log g = 7.60) and the highest N(R I) (T eff = 41,000 K and log g = 7.90). We find logiV(H I) = 
17.85 ± 0.03 ± 0.02, where the first uncertainties are statistical and the second are due to the 
different stellar models. We use a conservative approach to combine the two uncertainties 
and adopt logiV(HI) = 17.85 ± 0.05. Our measured iV(HI) is in good agreement with 
previous determinations using Lya data obtained with the IUE satellite (log iV(H I) = 17.92 
± 0.10, ?) and EUVE data (???, log iV(HI) = 17.90, log JV(HI) = 17.89 ± 0.04, and 
log iV(HI) = 17.90 ± 0.04, respectively). 



1213 1214 1215 1216 1217 1218 1219 
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Fig. 3. — Lya absorption for WD 1254+223. The histogram shows the coadded and merged 
STIS E140H observations of this line of sight; the thick solid line shows the best fit stellar 
model (T cS = 40,000 K and log g = 7.75). The other two models plotted produce the lowest 
(T eff = 39,000 K, log g = 7.60, dashed line) and highest (T cff = 41,000 K, log g = 7.90, 
dashed-dotted line) iV(HI). 
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4.3. WD 2152-548 

4.3.1. Stellar Model 

Determinations of T cS for this white dwarf range from T cS = 44,000 K (?) to T cff = 46,000 
K (?) while log g ranges from 7.8 to 7.9 (??). To take into account these different deter- 
minations of T c ff and log g as well as their associated uncertainties we constructed a grid 
of nine pure-hydrogen NLTE stellar models corresponding to T c g= (44,000;46,000;48,000) K 
and log g= (7.70;7.85;8.00). Figure 4 displays the stellar model with T c g = 46,000 K and 
log g = 7.85 overplotted on the STIS E140M data for the Lya region. 

4.3.2. Lya Profile Fitting 

For each of the nine stellar models computed we normalize the data by the stellar model 
prior to profile fitting the Lya HI transition (no DI is included in the profile fit). We use 
a single-Gaussian LSF with a FWHM of 1.3 pixels. In addition, a second-order polynomial 
is used to model the continuum after normalization by the model atmosphere, to take into 
account uncertainties in the instrument sensitivity. ? determined that «ph — ^ism ~ — 12 km 
s -1 . Profile fits taking into account the velocity shift between photospheric and interstellar 
absorptions have shown that this shift is negligible and thus does not affect the interstellar 
AT(HI) determination. 

As with WD 0050— 332, profile fitting of the Lya transition alone without constraining 
the b value leads to unreasonable estimates of iV and b. To determine the b value we use the 
FUSE data from the SiC IB channel which covers the Lyman edge. By fitting the H I Lyman 
series transitions between 913.6 A and 915.5 A with the continuum fixed, we determine &hi 
= 8.5 t. l\\ km s _1 . This b value is then used as a constraint in determining AT(HI) using 
the Lya transition. We find log iV(HI)= 18.84 t o!°4 ^ 0-01- The uncertainties quoted are 
statistical (associated with using the best fit model, T e g- = 46,000 K; log g = 7.85), and 
systematic, respectively. Using a conservative approach to combine the two uncertainties we 
adopt log iV(HI)= 18.84 + g;gf. Similarly to WD 0050-332 our derived b H1 implies T max is 
in disagreement with the mean temperature derived by ?. To determine the influence of &hi 
in the derived iV(HI) we performed multiple fits using the best fit stellar model and fixing 
6 H i at different values. We find that iV(HI) agrees with our adopted value at the la level for 
fr H i < 25 km s -1 . Hence our derived log iV(HI)= 18.84 + ^ is not sensitive to the assumed 
&hi, for reasonable values of 6hi- 

Our iV(H I) determined from Lya is in excellent agreement with that determined by ? 
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from EUVE data, who find log iV(HI) = 18.84 + oil, using pure-hydrogen LTE models with 
T eff = 44,300 K and log g = 7.91. 

5. The Abundance of Oxygen in the Local Bubble 

The first part of Table 6 presents the column densities and ratios derived in this work. 
In the second part of Table 6 we present a compilation of AT(OI) and iV(HI) (and their 
respective ratios) from the literature, for sightlines probing only the Local Bubble. Since 
the ionization balance for both O and H is linked by resonant charge exchange reactions 
(see ?, and references therein) and there is no H 2 along these sightlines, OI/HI should be 
representative of O/H along these lines of sight. Hence, in Table 6, O I/HI has been replaced 
by O/H. The last column in the table provides information about the sources of the column 
densities. The choice of the stars listed in the second part of Table 6 is based on the contours 
of the Local Bubble which have been mapped by ? and ? using absorption-line studies of 
Nal. They found a cavity with a radius between 65 and 250 pc, that is delineated by a sharp 
gradient in the neutral gas column density with increasing radius, a dense neutral gas "wall" . 
? quote an equivalent width of 20 mA of Nal A5891.59 as corresponding to log iV(HI) ~ 
19.3. This iV(HI) is used as a rough estimate of the Local Bubble boundary throughout this 
work and as such we only searched the literature for stars with log iV(HI) < 19.3. We note 
that in the literature there are other stars that satisfy this criterion (e.g., ?). However, for 
these stars, only the strongly saturated A1302 01 transition was used to determine iV(OI). 
We suspect that those column densities are not as reliable as the ones derived here (see 
note in Table 6) and as such we chose not to include them in Table 6. Hebrard et al. (in 
preparation) have also demonstrated, based on the constancy of the D/O ratio in the Local 
Bubble (?), that N(Ol) derived by ? are unreliable. 

? have also derived O/H along the sightline to WD 1029+537, which probes only gas 
in the Local Bubble. We chose not to include their value in Table 6 because we believe that 
their iV(OI) determination might be based on saturated lines (A1039 and A988.8/988.6). In 
addition, the inconsistency between the derived b- values for the O I transitions used in their 
analysis, 6 i = 11.2 ± 1 km s" 1 (A1039) and 6 OI = 30 ± 5 km s" 1 (A988.8/988.6), leads us 
to believe that N(0 1) is not reliable. 

All stars listed in Table 6 are white dwarfs, but Capella and a Vir. White dwarfs are 
particularly useful for studying the ISM in the Local Bubble, because they provide a source 
of light that covers a wide spectral window ranging from the extreme ultraviolet to the op- 
tical. The relatively simple spectra of white dwarfs are also another asset when analyzing 
absorption lines formed in the ISM. The intense gravitational field that characterizes white 
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Fig. 4.— Interstellar H I Lya toward WD 2152-548. The histogram shows the STIS E140M 
observations of this line of sight; the thick solid line shows the adopted stellar model, corre- 
sponding to T eS = 46,000 K and log g = 7.85. Also displayed are the models that produce 
the lowest (T eff = 44,000 K, log g = 7.70, dashed line) and highest (T eff = 48,000 K, log g = 
8.00, dashed-dotted line) iV(HI). Geocoronal emission is annotated with ©. 
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dwarfs leads to the downward diffusion of elements heavier than the dominant constituent of 
the atmosphere (either H or He), leaving the atmosphere devoid of heavy elements. There- 
fore, a typical spectrum of white dwarf shows a smooth continuum with only broad hydrogen 
or helium lines. The simplicity of white dwarf atmospheres minimizes the problem of blend 
between ISM and photospheric lines. Observations show, however, that traces of metals 
are detected in some white dwarf atmospheres. The presence of these heavy elements is at- 
tributed to radiative forces or weak stellar winds that counteract efficiently the gravitational 
settling. For stars showing the presence of heavy elements, a comparison between a stellar 
model and the observations allows us to identify ISM lines that are blended with photo- 
spheric lines. Consequently, no O I lines selected for determining the O I column densities 
presented here are blended with photospheric lines. 

The H I column densities along the lines of sight toward the white dwarfs listed in Table 6 
are obtained by fitting the damping wings of the interstellar H I Lya profile. The only excep- 
tion is WD 1634—573 for which iV(H I) has been obtained by fitting the Extreme Ultraviolet 
Explorer (EUVE) spectrum with a stellar model and by adjusting the HI, He I and Hell 
column densities (see ?, and references therein). Figures 2, 3, and 4 illustrate that the damp- 
ing wings of the ISM H I Lya profile span a small portion of the photospheric absorption, 
hence accurate stellar atmosphere models are needed to reproduce the stellar contribution 
to the Lya profile, in order to obtain the ISM profile. This is accomplished by computing 
stellar atmosphere models under the assumption of non-LTE and by incorporating detailed 
photospheric Lya profiles. We note that ? have suggested that for lines of sight with iV(H I) 
< 10 19 cm" 2 , HI column densities based on the Lya transition might be overestimated by 
as much as 20%, due to the presence of hot HI components. However, ? compared iV(HI) 
values derived from HST and EUVE data and found them to be in good agreement. 

The four ratios derived here are in good agreement with previously determined Local 
Bubble ratios, presented in the second part of Table 6. The O I transitions used to derive 
N(0 1) have been used in other works (see for example ?). In addition tests have shown that 
even though some O I /-values might be slightly inaccurate, the errors should decrease in 
the final result when multiple 01 transitions are used (?). HI column densities have been 
determined from the Lya transition numerous times (see for example the iV(H I) references 
in the second part of Table 6). The good agreement between our iV(HI) and the ones 
determined with a completely different method using EUVE data, gives us confidence in 
our iV(H I) determinations. 

Using all the ratios in Table 6 (except the one for which a lower limit is quoted) we 
derive the weighted average (0/H) L b = (3.45 ± 0.19)xlCT 4 (la in the mean; average of 
individual uncertainties used to compute a) with xl — 11 f° r 10 degrees of freedom. The 
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dispersion around the mean as measured by the square root of the weighted average variance 
is 0.65 xlO~ 4 . Figure 5 displays the ratios used to compute the weighted average above, as 
well as the derived (0/H) LB (dashed line) and the Solar value from ? (dotted line). 

To determine whether any of the ratios dominates the x 2 statistics we determined the 
X 2 , repeatedly, with all the ratios but one. We found that the ratio that most influences the 
xl is the one derived towards WD 1254+223, O/H = (2.51 1 ^xlO -4 . When this ratio is 
excluded from the weighted average we find O/H = (3.66 ± 0.21) xl0~ 4 with xl — 0.6 for v 
= 9. The corresponding square root of the weighted average variance is 0.49 x 10~ 4 . However, 
this new ratio is still consistent at the la level with the weighted average determined with 
the first eleven ratios in Table 6, O/H = (3.45 ± 0.19)xl0 -4 . Because the uncertainties vary 
from sightline to sightline, some of the measurements have a small influence on the value of 
the weighted mean, the uncertainty in the mean and the xl- Thus, the sightlines with the 
largest uncertainties, G191— B2B, a Vir, and Capella, change the weighted mean by only 
~2% and decrease the uncertainty in the mean by ~3%. 

Based on our new and previous measurements of O/H, we conclude that this ratio is 
constant in the Local Bubble, and that (0/H) LB = (3.45 ± 0.19) xlO" 4 . 

6. Discussion 

Our (0/H) LB = (3.45 ± 0.19) xlO" 4 is virtually identical to O/H = (3.43 ± 0.15) xl0~ 4 
(?, updated /-value) derived using a single transition (A1356), for 13 lines of sight probing 
gas within 1500 pc (most within 500 pc), with 20.18 < log JV(HI) < 21.28. The /-value 
for A1356 (/ = 1.16xl0~ 6 ) as suggested by ? and ? is based on the mean of the lifetimes 
measured by several workers (????) and the theoretical branching ratio determined by ?. 
Our result is also in agreement with O/H = (3.47 ± 0.07) xl0~ 4 derived by ? for sightlines 
within 800 pc with average density n H < 1.0 cm -3 . 

The /-values of the O I transitions used to determine all the iV(0 1) summarized in Table 
6 have not been experimentally determined and as such are subject to larger uncertainties. 
However, the similarity of the O/H ratio determined here with that of ? seems to suggest 
that the net effect of the uncertain /-values on iV(OI) is not significant (but see ?, for two 
/-values, not used here, that are clearly overestimated). 

By comparing our gas-phase ISM oxygen abundance with the Solar value, and by as- 
suming that the Solar value represents the total (gas + dust) abundance of oxygen, we can 
estimate the amount of oxygen that can be incorporated into dust. The commonly adopted 
Solar oxygen abundance compiled by ?, (O/H) ~8.5xl0~ 4 , has been revised downwards 
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Fig. 5.— O/H ratios used in deriving (0/H) LB = (3.45 ± 0.19) xKT 4 (dashed line), as a 
function of the hydrogen column density. The Solar value (O/H)© = 4.6 xlO -4 from ? is 
indicated by a dotted line. 
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during the last few years. ? determined (O/H) ~5.5xl0" 4 and ? (O/H) ~4.9xlO" 4 . 
More recently ? found (O/H) = (4.60 ± 0.50) xlO" 4 and ? determined (O/H) ~ 4xl0~ 4 . 
This downward trend in the Solar oxygen abundance implies that a smaller amount of oxy- 
gen can be incorporated into dust. If we adopt the new solar oxygen abundance of ? as 
representative of the Solar abundance of oxygen, then ~25% of the oxygen in the Local 
Bubble could be depleted onto dust grains. However, this should only be a lower limit to 
the amount of oxygen depleted in grains since stellar nucleosynthesis is supposed to increase 
the oxygen abundance from the protosolar value. 

In contrast to the D/O and D/H ratios which seem to have different values outside and 
inside of the Local Bubble (??), our results are consistent with a constant and single- valued 
O/H ratio, both inside and outside the Local Bubble, up to log iV(HI) ~ 21. 

The composition of the Local Interstellar Cloud (LIC) has been studied using pickup 
ions observed inside the heliosphere (?). These authors report atomic densities that yield 
OI/HI = 4.32 xlO -4 compared to our value of 3.45 xlO -4 . ? have shown in a study of D 
I/O I ratios for 14 sightlines that the variation in the LB was small, consistent with the 
measurement uncertainties. Unless one assumes that O I and D I abundances are correlated, 
contrary to the usual assumptions of astration, both the O I and D I abundances do not vary 
by large amounts. Thus, the heliospheric results are indicative of abundances in the LB. 
The uncertainties in the densities at the termination shock and filtration factors reported 
by ? imply that the uncertainty in the ratio above is 21%. The agreement within the 
uncertainties between the heliospheric value and the spectroscopic value, which were obtained 
by very different techniques, is reassuring and indicates that large systematic errors (such as 
incorrect oscillator strengths) are not significant. 



6.1. Implications of (0/H) LB on the Galactic Infall Scenario 

The similarity of the value of O/H in the Local Bubble reported here with the values 
measured at larger distances (???) is striking and may provide a limit on the infall of low 
metallicity gas into this part of the Milky Way. The quantity of gas above and below the 
Milky Way plane is likely to be much less in the vicinity of the LB compared to that in 
the LISM studied by ?. The average densities, defined as the column density divided by 
distance, for the LB sightlines reported here ranges from 0.003 to 0.08 cm -3 with a median 
value of 0.02 cm -3 (mean value is 0.03 cm -3 ). In comparison, the sightlines used by ? to 
obtain their average value, for d < 800 pc and riu < 1.0 cm -3 , range from 0.09 to 0.91 cm -3 
with a median value of 0.26 cm -3 (mean is 0.34 cm -3 ). Although the average density is not 
a direct measure of the gas in a column height above or below the plane, it is a qualitative 
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indicator of the amount of gas in each region. Thus, for similar infall rates, dilution of metals 
will be larger in the LB compared to the LISM. Using the uncertainties in the mean in this 
work and that of ? as limits, the maximum possible difference between the LB and the 
LISM is 0.32 xl(T 4 , 9.3% of the LB value of O I/H I. If the quantity of gas in the vicinity of 
the LB is a factor of 3 below that beyond the LB, then by using simple dilution arguments 
and assuming infalling material with properties similar to that studied by ? in Complex 
C with OI/HI, 0.17 that of solar ((O/H) = 4.6 xl0~ 4 from ?), we find the fraction of 
infalling gas within a mixing time to be no more than about 15%. This modest amount 
of infall would have a small effect on the D I abundance in the LB. Using the value of D 
I/H I = 2.2 xlO" 5 reported for the Complex C material by ?, we find that if the present 
LB value is 1.56 xlO -5 (?), the D/H ratio would have been increased by infall from a value 
of 1.46 xl0~ 5 , a negligible difference. It is unlikely that such modest infalls can be used to 
explain the difference between the LB D/O and D/H values and the lower values found for 
more distant sightlines with typically d > 500 pc and log iV(HI) > 20.5 (???). If we had 
compared our value with that of ?, which has smaller uncertainties, the infall limit would 
be even smaller. Finally, we note that a discussion of infall without including mixing and 
chemical evolution is extremely simplistic and is included here only to indicate the magnitude 
of the effects implied by the small diference between the LB and the LISM values. 



6.2. Implications of (0/H) LB on the Population of O-bearing Grains 

The similarity between the O/H values for the LB and LISM also allows us to speculate 
about the nature of the grain population in the two environments. The sightlines studied 
here, together with the ones from ? and the ones from ? discussed above, probe different 
average line of sight densities, n H . As discussed above, the average densities for the Local 
Bubble sightlines presented here are all <0.1 cm -3 , whereas the densities beyond the LB 
are higher. The constancy of the O/H ratio across this density range may indicate that the 
population of oxygen bearing grains is similar in the different environments. In contrast, 
many other elements vary with density (see the review by ?). ? discussed how the iron 
depletion varies between sightlines with different average gas density, using a simple model 
proposed by ?. For densities less than about 0.1 cm~ 3 , they pointed out that the sightline 
would only pass through clouds (or density fluctuations) corresponding to warm diffuse gas. 
With increasing density to a few tenths of atoms per cm 3 , sightlines would also pass through 
cool diffuse clouds, and at even higher densities, cool dense clouds with much lower spatial 
frequencies, would be sampled. Fitting their data to a model, they found that Fe is a factor 
of four less depleted in warm diffuse clouds than in cool diffuse clouds, and suggested that 
the dust distribution was likely different in the two different diffuse environments. Even 
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though Fe/H varies in this «h range, O/H does not. It is possible that this change in the Fe 
depletion is due to the presence of Fe in grain mantles that are easily created and destroyed 
in the interstellar medium. It is also possible that the oxygen bearing grains and the grains 
responsible for the changes in the iron depletion are different populations. 

A related question may be why does O/H remain invariant beyond the LB, while D/O 
and D/H do not (??, and references therein). In a recent examination of the question whether 
dust can explain the measured variations in the D/H ratio, ? points out that outside dark 
clouds, the oxygen in grains is thought to be mostly in the form of Mg :r Fe2_ x Si04. On this ba- 
sis and using the interstellar abundance of Si (?), he estimated (0/H) gas « 0.76x(O/H) to tai, 
a value similar to that found here and by others (??). If the O-bearing grains are as described 
above and are robust against destruction in the warm diffuse and cold diffuse media, what 
are the Fe-bearing grains, responsible for the variable depletion? These grains may be the 
aromatic carbonaceous grains ? suggests are responsible for depleting deuterium and hence 
the variation in D/O and D/H ratios (?, and references therein). ? have already suggested 
that ~60% of cosmic Fe could be associated with a population of small carbonaceous grains. 

? suggests that if the scenario discussed above is correct, the destruction of these grains 
by shocks with v s > 100 km s _1 will release both D and Fe into the interstellar medium. J. J. 
Linsky (2005, in preparation) shows that the D and Fe abundances appear to be correlated, 
and C. Oliveira (2005, in preparation) has shown that D/O and D/H decrease with increasing 
density, for ran > 0.1 cm~ 3 . Thus it is possible that the reason why O does not vary, unlike 
Fe and possibly D, is the robust nature of the O-bearing grains, as described above. If 
future studies strengthen the evidence for dust depletion of D, they will also strengthen the 
hypothesis of distinct grain populations. 

7. Summary 

We have used FUSE data in conjunction with STIS observations and previously pub- 
lished OI and HI column densities to determine the O/H ratio in the Local Bubble. We 
find the weighted average (0/H)lb = (3.45 ± 0.19) xlO -4 (la in the mean) for a total of 
eleven sightlines. Our result is in agreement with previous determinations of this ratio along 
sightlines probing gas well outside the Local Bubble (??). 

Based on the similarity of the two ratios above, our results suggest that the net effect of 
the uncertain FU SE bandpass O I /-values is not significant in the determination of N(0 1) 
for these LB sightlines. 

Taking into account the latest determinations of the Solar abundance of oxygen (??), 
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our results indicate that ~25% of the oxygen in the Local Bubble is tied up in dust grains. 

In conjunction with other studies (??) our work suggests limits on infall of low metal- 
licity gas in this part of the galaxy. It also suggests that the populations of O-bearing grains 
are similar in both warm and cold diffuse clouds. 

This work is based on data obtained for the Guaranteed Time Team by the NASA- 
CNES-CSA FUSE mission operated by The Johns Hopkins University. Financial support to 
U. S. participants has been provided in part by NASA contract NAS5-32985 to Johns Hopkins 
University. Based on observations made with the NASA/ESA Hubble Space Telescope, 
obtained from the Data Archive at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc. under NASA contract 
NAS 5-26555. The profile fitting procedure, Owens. f, used in this work was developed by 
M. Lemoine and the French FUSE Team. We thank Nicolas Lehner for helpful discussions 
regarding iV(OI). 
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Table 1. Stellar Properties 



Star 




T e ff 


log g 


d 


I 


b 






(K) 


(cm" 2 ) 


(pc) 






WD 0050- 


-332 


36320 ± 200 


7.87 ± 0.04 


58 


299.14 


-84.12 


WD 0232H 


-035 a 


55,000-57,000 


7.3-7.5 


70 


165.97 


-50.27 


WD 1254H 


-223 


39600 ± 200 


7.80 ± 0.03 


67 


317.26 


+84.75 


WD 2152- 


-548 


45800 ± 600 


7.78 ± 0.07 


128 


339.73 


-48.06 



Note. - T e fj, log g and d are from ? unless noted otherwise. I and 
b are from the SIMBAD database. All white dwarfs have spectral type 
DA. 

a T ef f, log g, and d from ?. 
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Table 2. Log of FUSE observations 



Star 


Program ID 


Aperture a 


Mode 


Time (ks) 


Date 


WD 0050-332 


M1010101 


L 


TTAG 


16.4 


2000 Jul 4 




P2042001 


L 


TTAG 


8.5 


2000 Dec 11 


WD 0232+035 


P1040501 


M 


HIST 


2.3 


2003 Dec 7 




P1040503 


M 


HIST 


16.1 


2004 Jan 2 




P1040504 


M 


HIST 


12.0 


2004 Jan 6 


WD 1254+223 


M1010401 


L 


TTAG 


6.3 


2000 Mar 6 




M1010402 


L 


TTAG 


12.2 


2000 Apr 29 




P2041801 


L 


TTAG 


9.9 


2001 Jan 28 




M1010403 


L 


TTAG 


8.1 


2001 Feb 7 


WD 2152-548 


P2042301 


M 


TTAG 


15.5 


2001 May 22 




M1051501 


L 


TTAG 


4.7 


2002 Sep 24 



a M and L stand for MDRS and LWRS, respectively. 



Table 3. Log of STIS observations 



Star 


Program ID 


Grating 


Exp. Time (ks) 


Date 


AA (A) 


WD 0050-332 


O4G101010 


E140H 


2.018 


1999 Jan 14 


1170-1372 




O4G101020 


E140H 


2.116 


1999 Jan 14 


1170-1372 


WD 1254+223 


O67K03010 


E140H 


1.851 


2001 March 14 


1140-1335 




O67K03020 


E140H 


2.632 


2001 March 14 


1140-1335 




O67K03030 


E140H 


2.632 


2001 March 14 


1140-1335 




O67K03040 


E140H 


1.950 


2001 March 14 


1140-1335 




O67K03050 


E140H 


2.632 


2001 March 14 


1140-1335 




O67K04010 


E140H 


1.851 


2001 March 22 


1140-1335 


WD 2152-548 


O4G105010 


E140M 


2.171 


1999 May 23 


1140-1735 



Note. — All the observations were performed through the 0.2" x 0.09" aperture. 
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Table 4. Equivalent Width Measurements (mA) 



A (A) 


/-value 


WD 0050-332 


WD 0232+035 


WD 1254+223 


WD 2152-548 


ioUz.iooO 


4.80xl0~ 2 


100.2 ± 1.9 


110.4 ± Z.U 


00. z ± o.U 


1 no A I r A 


1039.2304 


9.07xl0~ 3 


49.1 ± 1.3 


48.2 ± 1.0 


13.7 ± 0.9 


63.1 ± 2.5 


976.4481 


3.31xl0~ 3 


25.5 ± 2.4 


25.4 ± 1.8 




39.2 ± 6.6 


971.738 a 


1.38xl0~ 2 


41.5 ± 5.4 


48.2 ± 2.8 




76.6 ± 11.2 


948.6855 


6.31xl0~ 3 


38.5 ± 4.3 


32.4 ± 2.5 




61.5 ± 10.9 


936.6295 


3.65xl0~ 3 


25.2 ± 3.9 


24.5 ± 1.9 




49.6 ± 6.6 


929.5168 


2.29xl0~ 3 


23.9 ± 4.3 


10.3 ± 1.5 




40.6 ± 8.6 


924.950 


1.54xl0~ 3 


19.3 ± 4.1 








Note. - 


Wavelengths and /-values are 


from ?. 







Triplet structure used, see reference above. 



Table 5. 01 Column Densities (Log) 



Method 


WD 0050-332 


WD 0232+035 


WD 1254+223 


WD 2152-548 


AOD 
COG 
PF 


15.15 t (924) 
14.98 ± 0.04 
15.09 ± 0.09 


15.01 ± 0.03 (976) 

14 Qfi + °- 04 

14 04 + o-os 
14. y4 _ 04 


14.22 ± 0.03 (1039) 
14.27 ± 0.04 


15.35 t Ho (936) 

1 r r 2 + 0.07 
IO. OZ _ 0.O6 

15 30 + a13 

IO. OU _ 0.09 


Adopted 


15.07 ± 0.09 


15.00 ± 0.06 


14 9c + 0.06 
14. ZO _ o.05 


15.45 ± 0.10 



Note. — AOD, COG, and PF stand for apparent optical depth, curve of growth, and profile 
fitting, respectively. The transition used to determine A^(OI) with AOD is indicated in (). 
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Table 6. Column Densities and O/H Ratios 



Star Distance (pc) logiV(HI) logiV(OI) O/HxlO 4 Ref(HI), Ref(OI) 



WD 0050-332 


58 


18.57 + g|« 


15.07 ± 0.09 


3 16 + a76 

°- iu - 0.65 


1, 1 


WD 0232+035 


70 


18.48 t gig* 


15.00 ± 0.06 


q qq + 0.59 
000 - 0.56 


2, 1 


WD 1254+223 


67 


17.85 ± 0.05 


14 OK + 0.06 


on + 0.48 
Z "< J± - 0.39 


1, 1 


WD 2152-548 


128 


18.84 ± 0.05 


15.45 ± 0.10 


4 07 + 116 

^• u ' - 0.95 


1, 1 


WD 0642-166 


2.6 


17.81 ± 0.11 


14.29 ± 0.05 


q no + 0.95 
°- uz - 0.75 


3,4 


HD 34029 


12.5 


18.22 ± 0.08 


15.02 ± 


6 31 + 3 - 69 

- 1.95 


5, 6 


WD 1634-573 


37 


18.85 ± 


15.51 ± 0.03 


4.57 ± 0.75 


7,8 


WD 1314+293 


68 


17.93 ± 0.03 


14.49 ± 0.04 


q fjq + 0.44 
O.UO _ 040 


9, 9 


WD 0501+527 


69 


18.18 ± 0.09 


14.86 ± 0.04 


A7Q + L19 
^- ' y - 0.99 


10, 10 


WD 2309+105 


79 


19.11 ± 0.03 


1 5 67 + °- 04 

10.0 1 _ o.o3 


q co + 0.44 
°- uo - 0.34 


11, 11 


HD 116658 


80 


19.00 ± 0.10 


15.58 ± 0.10 


3.80 + If, 


12, 12 


WD 1029+537 a 


132 


18.62 ± 0.05 


>14.93 


>2.04 


13, 1 



a iV(OI) from ? determined from the A1302 transition is in disagreement with the lower limit 
determined from the weaker A1039 transition observed by FUSE which we report here. 



References. — (1) This work; (2) ?; (3) ?; (4) ?; (5) ?; (6) ?; (7) ?; (8) ?; (9) ?; (10) ?; (11) ?; 
(12) ?; (13) ?. 



